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ABSTRACT 


Magnetars!? 


are young, highly magnetized neutron stars that produce extremely rare giant 
flares of gamma-rays, the most luminous astrophysical phenomena in our Galaxy. The detection of 


these flares from outside the Local Group of galaxies has been predicted ^?, with just two candidates 


so far°-!°. Here we report on the extremely bright gamma-ray flare GRB 200415A of April 15, 
2020, which we localize, using the Interplanetary Network, to a tiny (20 sq. arcmin) area on the 
celestial sphere, that overlaps the central region of the Sculptor galaxy at ~ 3.5 Mpc from the 
Milky Way. From the Konus-Wind detections, we find a striking similarity between GRB 200415A 
and GRB 051103, the even more energetic flare that presumably originated from the M81/M82 
group of galaxies at nearly the same distance (3.6 Mpc). Both bursts display a sharp, millisecond- 
scale, hard-spectrum initial pulse, followed by an approximately 0.2 s long steadily fading and 
softening tail. Apart from the huge initial pulses of magnetar giant flares, no astrophysical signal 
with this combination of temporal and spectral properties and implied energy has been reported 
previously. At the inferred distances, the energy released in both flares is on par with that of the 
December 27, 2004 superflare*!!? from the Galactic magnetar SGR 1806-20, but with a higher 
peak luminosity. Taken all together, this makes GRB 200415A and its twin GRB 051103 the most 
significant candidates for extragalactic magnetar giant flares, both a factor of = 5 more luminous 
than the brightest Galactic magnetar flare observed previously, thus providing an important step 


towards a better understanding of this fascinating phenomenon. 


MAIN TEXT 


Magnetars?? are a special rare class of neutron stars with strong magnetic fields! (B ~ 10!^ — 10? G). 
Some magnetars (or Soft-Gamma Repeaters, SGRs) exhibit bursting emission in hard X-rays/soft gamma- 
rays. During the active stage, which may last from several days to years, SGRs emit randomly occurring 
short (from milliseconds to seconds long) hard X-ray bursts with peak luminosities of ~ 1038-10% erg s7!. 
Much more rarely, perhaps several times during the SGR stage!!, assumed to last up to ~ 10? years which 
corresponds to the SGR magnetic field decay time-scale, a magnetar may emit a giant flare (GF) with the 
sudden release of an enormous amount of energy in the form of gamma-rays ~ (0.01 — 1) x 1076 erg. A GF 
displays a short (fraction of a second) initial pulse of radiation with a sharp rise and a more shallow decay 
which evolves into a soft spectrum, long-duration decaying tail modulated with the neutron star rotation 


period. Only a dozen burst-emitting magnetars in our Galaxy and Large Magellanic Cloud are known so 


far and only three of them produced a GF during the active stage!*. Due to the enormous luminosity of the 


initial pulse, GFs can be detected from magnetars in nearby galaxies up to tens of megaparsecs away*?. In 
this case, the initial pulse may mimic a short gamma-ray burst (GRB) produced by the merger of a binary 
compact object (two neutron stars or a neutron star and a black hole) in galaxies at cosmological distances. 
Neither the rotationally modulated tail, nor the short weaker X-ray bursts preceding and following a GF 
can be observed even from nearby galaxies by the current wide-field gamma-ray monitors. Thus, the main 
evidence supporting the magnetar nature of such bursts are spatial coincidence with a galaxy, consistency 
of the burst light curve and energetics with known GFs, and non-detection of a gravitational wave signal 
indicative of a short GRB. Up to now only two short GRBs were proposed to be magnetar giant flares 
outside our Galaxy and its satellites (extragalactic magnetar GFs, eMGFs): GRB 051103, associated 
with the M81/M82 group of galaxies®® at a distance Dus /M82 = 3.6 Mpc!* and GRB 070201, in the 
Andromeda galaxy (M31)?!° at Dy31 = 0.77 Mpc!. Detailed searches for GFs from nearby galaxies 


16,17 


have not revealed any other credible candidates and put an upper limit on the fraction of GFs in the 


short GRB population less than ~ 10%. 

On 15 April 2020 the extremely bright, short GRB 200415A occurred at 08:48:06 UT at the Earth, and 
was detected by five missions of the Interplanetary network of gamma-ray detectors (IPN, see Methods). 
We report here the final ~ 20 arcmin? burst localisation by the IPN (see Methods), which overlaps the 
central part of the nearby Sculptor galaxy, also designated as NGC 253, at a distance DyGc253 = 3.5 Mpe!’ 
(Figure 1). The chance occurrence for GRB 200415A to be spatially consistent with a nearby galaxy likely 
to produce detectable eMGFs is approximately 1 in 2000001”. 

A preliminary analysis of the Konus-Wind (KW) detection? of GRB 200415A revealed a remarkable 
similarity between this burst and GRB 051103, historically the first eMGF candidate. To explore this 
similarity further, we performed a detailed comparative study of their temporal and spectral properties 
(Methods); the results are in general agreement with those previously reported on the latter event^?. 

As observed by KW, the 2 ms light curves of GRB 200415A and GRB 051103 start with the fast, 
(~ 2 ms) rise of a narrow (~ 4 ms) initial pulse (IP), which is followed by an exponentially decaying phase 
with count-rate e-folding time Ter ~ 50 ms (Figure 2). The total burst durations are 0.138 s (GRB 200415A) 
and 0.324 s (GRB 051103), and the values of Tan (the duration of the time interval which contains the 


central 90% of the total burst count fluence) are 0.100 X 0.014 s and 0.138 +0.005 s, respectively (hereafter 
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all the quoted uncertainties are at the 68% confidence level). Although the peak count rates, reached in the 
first 2 ms of the IPs, are very similar, ~ (1.5-1.7) x 10? s-!, the photon flux over the entire extent of the 
decaying phase is about twice as high in GRB 051103 as in GRB 200415A. 

The burst hardness (the ratio between the 390-1600 keV and 90-390 keV count rates) rapidly 
increases during the initial pulses, peaks during the following ~ 8 ms, and then gradually decays with 
the burst count rate. Our spectral analysis (Methods) shows that starting from the rise of the IP and up 
to ~ To + 100 ms, the energy spectra of both bursts are well described by a cut-off power law function 
(CPL; ~ ET exp(—E(a -- 2)/ Ey)). The temporal evolution of the spectra is illustrated in Figure 2 which 
shows the behaviour of the CPL model parameters: the peak energy of the EF (E) spectrum Ep (panel c) 
and photon power-law index G (panel d). The initial narrow pulses of both bursts are characterized by 
Ep ~ 1.2 MeV, with a ~ —0.6 for GRB 200415A and a much harder œ ~ —0.1 for GRB 051103. This Ep 
was the highest reached in the entire event for GRB 200415A, while the hardest emission in GRB 051103 
(Ep ~ 3 MeV, G ~ 0.2) was reached during the subsequent ~ 30 ms. A non-thermal CPL model,with Ep 
decaying nearly exponentially, adequately describes burst spectra up to ~ 7o + 100 ms. Afterwards, the 
very hard photon index @ becomes poorly constrained and, simultaneously, the emission spectrum can 
be described by a blackbody function (with a temperature kT ~ 70—100 keV), which is excluded by our 
analysis at the initial stages of the bursts. 

Panel b of Figure 2 shows the temporal evolution of the 20 keV-10 MeV energy flux. For both bursts, 
the flux peaks in the initial spike and, starting from ~ Tọ + 50 ms, decays with tg ~ 30 ms. In accordance 
with the burst similarities in peak count rate and energy spectrum, measured in the IPs, their 4 ms peak flux 
estimates also agree within errors: 0.96772 x 1073 erg cm? s^! and 1.15*05; x 1073 erg cm? s^! for 
GRB 200415A and GRB 051103, respectively. 

The time-integrated spectra of both flares, measured from To to 7o + 0.192 s, are best described by a 
sum of non-thermal (CPL) and thermal (BB) components. The burst 20 keV-10 MeV fluences are 85m 
(34.3*50) x 107 erg cm? for GRB 200415A (GRB 051103), with blackbody component contributions 
of ~ 14% and ~ 9%, respectively. The contribution of the initial short spike to the total burst fluence is 


about 45% for GRB 2004154, and just 13% for GRB 051103. 


Thus, the extremely bright short GRB 2004154, which strong evidence suggests is associated with 
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the Sculptor galaxy, is strikingly similar to GRB 051103, that presumably originated from the M81/M82 
group of galaxies at nearly the same distance, in terms of light curve morphology, spectral behavior, and 
observed peak energy flux. A lightcurve with a bright, millisecond scale initial pulse followed by an 
exponentially decaying emission is quite unusual for short cosmological GRBs; none of = 500 short 
bursts detected by Konus- Wind in more than 25 years of observations displays such a shape^?!. On the 
other hand, this pattern was observed in two Galactic magnetar giant flares, from SGR 1900+14°%?3 and 
SGR 1806-20??-?*, Furthermore, higher time resolution light curves of GRB 200415A from Swift-BAT 
and Fermi-GBM~ show an initial short (< 1 ms) sub-peak, followed by a sharp decrease for ~ 1 ms, 
before the main part of the peak. This pattern is also seen in the SGR 1806-20 giant flare!! and may 
be a general property of magnetar GFs that can be used to identify them within the short GRB sample. 
Thus, the interpretation of both GRB 200415A and GRB 051103 as magnetar giant flares is strongly 
suggested, with additional support from the non-detection of an accompanying gravitational wave signal 


for GRB 05110370 (there is no sensitive coverage by a gravitational wave detector for GRB 200415A). 


At source distances Dycc253 = 3.5 Mpc and Dmg1 = 3.6 Mpc, the characteristic radius of the emission 
region, estimated from the blackbody spectral fits, is R ~ 20 — 40 km, of the same order as the radius of a 
neutron star or its magnetosphere. The implied isotropic-equivalent energy release in y-rays, Eiso, is ~ 1.3 
(~ 5.3) x10* erg, and the isotropic-equivalent peak luminosity, Liso, is ~ 1.4 (~ 1.8) x 10*8 erg s^! for 
GRB 200415A (GRB 051103). Thus, the total energies released in both flares are comparable with that 
estimated for the most energetic flare from a Galactic magnetar ~ 2.3 x 1049 erg’, but at much higher peak 
luminosity. Taken all together, this makes GRB 200415A and its twin GRB 051103 the most significant 
candidates for extragalactic magnetar giant flares, both a factor of = 5 more luminous than any Galactic 
magnetar flare observed previously’. Assuming the same spectra and energetics, similar events can be 
detected with KW from distances up to ~16 Mpc. The properties of GRB 200415A and GRB 051103 are 
summarized in Extended Data Table 1. 

Despite the strong evidence in favor of the magnetar GF nature of GRB 200415A and GRB 051103, it 
cannot completely be ruled out that they might belong to an as-yet undiscovered branch of the cosmological 
short GRB population. For the observed energy fluence, and assuming a cosmological redshift z of 0.05-1, 


GRB 200415A is consistent with the Konus-Wind sample of short GRBs with known redshifts?^?5 in 
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terms of a hardness-intensity relation in the cosmological rest frame (see Extended Data Figure 1). In the 
case of GRB 051103, the implied short GRB redshift is z ~ 1, and intrinsic Ep ~ 5 MeV. 

GRB 200415A, detected in the wide energy range ~ 10 keV up to GeV, by multiple space-based 
observatories, is to date the most comprehensively studied initial phase of MGF (the pulsating tail is 
beyond the sensitivity of current instruments from 3.5 Mpc’). The detection of extragalactic magnetar 
giant flares facilitates the study of emission processes on millisecond and sub-millisecond time scales? 
inaccessible to galactic events, which saturate almost all gamma-ray detectors. On the time scale accessible 
to KW (2 2 ms), the single-peaked GRB 200425A and GRB 051103 are clearly different from the third 
known eMGF candidate, GRB 070201°, with a highly variable emission during the first ~ 50 ms. This 
suggests that the physical processes behind the emission in MGF initial pulses may develop on temporal 
scales spanning more than an order of magnitude. The exponential decay in both the flux and Æp is also 
observed in the GBM/BAT with higher statistical quality, which is used to interpret the physical emission 
process and environment of the likely neutron star progenitor”. The recognition of GRB 200425A as 
an eMGF in the Sculptor galaxy has provided a clue for the interpretation of the delayed GeV photons 
detected by Fermi-LAT as the emission from interaction of a magnetar ultra-relativistic outflow with 
environmental gas??. Finally, the detection and the subsequent retrospective eMGF search!? provide an 


excellent proxy for the magnetar population in nearby galaxies. 
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Figure 1. The final IPN localization of GRB 200415A superimposed on an image of the Sculptor galaxy 
from the GALEX survey (1750-2800 À; see Methods). The localization is defined by the 4.73 arcmin 
wide Wind-Odyssey and 3.58 arcmin wide Wind-Fermi annuli. The IPN error box (purple parallelogram) 
is shown along with the 20 arcmin? 30 error ellipse (shown in magenta) for the position. The coordinates 
are J2000. 
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Figure 2. Time histories of the twin GRB 051103 (red) and GRB 200415A (blue) and evolution of their 
spectral parameters. All times are given relative to the KW trigger time 7p. Panel a shows the burst time 
histories as recorded by Konus-Wind. Both events start with a sharp rise of an exceptionally bright, narrow 
(4 ms) initial pulse (IP), whose time interval is indicated by a green shaded area in each panel, followed by 
an exponential decay with Ter ~ 50 ms (dotted lines). Panel b shows the evolution of the 20 keV-10 MeV 
energy flux, which, in both cases, peaks in the IP and, starting from ~ Tọ + 50 ms, decays with 

taux ^v 30 ms. The temporal evolution of the emission spectra is illustrated with the behavior of CPL 
model best-fit parameters: peak energy Ep (panel c) and photon power-law index & (panel d). Both bursts 
are characterized by Ep ~ 1.2 MeV in the IP, which is the hardest part of GRB 2004154, while the 
hardest emission in GRB 051103 (with Ep ~ 3 MeV) was detected during the subsequent ~ 30 ms. A 
non-thermal CPL model adequately describes burst spectra up to ^ 7o + 100 ms; afterwards, the hard 
power-law photon index œ becomes poorly constrained and, simultaneously, the emission spectrum can be 
described by a blackbody function with kT ~ 70—100 keV. Vertical error bars indicate 68% confidence 
intervals, and horizontal error bars indicate the duration of the interval. 
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Methods 


Observations 


GRB 200415A occurred at 08:48:06 UT at the Earth, and was detected by the Gamma-ray burst monitor 
(GBM) on-board the Fermi Gamma-ray Space Telescope Mission?? and the Burst Alert Telescope (BAT, 
outside the coded field of view) on-board the Neil Gehrels Swift Observatory?! in low Earth orbits; the 
SPI telescope anticoincidence system** (SPI-ACS) and Pixellated Imaging Caesium Iodide Telescope??? 
(IBIS-PICSIT) instruments on-board the International Gamma-Ray Astrophysics Laboratory (INTEGRAL) 
in an eccentric Earth orbit at 0.44 light seconds from Earth; the Mars-Odyssey High Energy Neutron 
Detector??? (HEND) in orbit around Mars at 672 It-s from Earth; the Konus- Wind?! (KW) gamma-ray 


38 in a Lissajous orbit at the L; libration point of the Sun-Earth system 


spectrometer on the Wind mission 
at a distance of 4.8 It-s; and the Modular X- and Gamma-Ray Sensor (MXGS) of The Atmosphere-Space 


Interactions Monitor ASIM?? on-board the International Space Station (not a part of the IPN). 


Two independent and consistent Fermi-GBM localizations^^^! RoboBA* and BALROG ^*^ were 
announced in GCN Circulars at 08:58 UT and 09:11 UT, respectively, each with 30 probability region 
covering more than 1000 deg?. Exploiting the difference in the arrival time of the gamma-ray signals at the 
four Interplanetary Network (IPN, see the next section) instruments (GBM, BAT, HEND, and SPI-ACS), 
a preliminary 1.5 deg? IPN error box was announced at 16:51 UT*, where it was pointed out that this 
event might be a magnetar GF in the Sculptor galaxy (NGC 253) located at Dygc253 = 3.5 Mpc. The 


Fermi-LAT localization, announced at 19:30 UT^Ó was consistent with the box. 


As soon as Konus-Wind data arrived the improved 274 arcmin” box (~ 20 times smaller than the initial 
box) was published on 16 April 16:16 UT“ which strengthened the association of the burst with the galaxy. 
The box was close to the 68% confidence contours of the updated Fermi-LAT localization*®, published on 
16 April 20:48 UT. The IPN localization was within about 37 deg of the Sun, making X-ray and optical 
follow-up observations challenging. The only optical observation of the Sculptor galaxy, resulting in an 


upper limit, was reported by the MASTER telescope network on 17 April^?. 
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Interplanetary network The Interplanetary network (IPN; http://www.ssl.berkeley.edu/ 
ipn3/) is a group of spacecraft orbiting the Earth and Mars equipped with gamma-ray burst detectors used 
to localize gamma-ray bursts??. When a GRB arrives at two spacecraft, it may be localized (triangulated) 
to an annulus on the sky, determined by the measured propagation time delay and spacecraft positions. 
Three spacecraft produce two possible locations (IPN error boxes). The ambiguity can be eliminated by 
the addition of a fourth, non-coplanar spacecraft, by the anisotropic response of KW, or by the GBM 
localization?! . 

The propagation time delay and its uncertainty are calculated by cross-correlation??. The systematic 
uncertainties are estimated using the comparison of IPN triangulations with precise GRB positions?!. 

The declared on-board clock accuracy of the spacecraft are: down to 1 us for Fermi; ~ 200 us for 
Swift; < 1 ms for Wind; ~ 100 us for INTEGRAL; for Mars-Odyssey an overall 30 systematic uncertainty 
which includes timing and other effects derived from IPN observations of precisely localized GRBs is 
better than 360 ms. The Wind clock drift information is provided at £tps: //pwgdata.sci.gsfc. 
nasa.gov/pub/wind_clock/. 

Near-Earth spacecraft ephemerides are derived from two-line elements (TLE) available at https: / / 
www. space-track.org using SGP8 model. The Wind predicted ephemerides data and their descrip- 


tion are available athttps://spdf.gsfc.nasa.gov/pub/data/wind/orbit/pre_or/ and 


https://cdaweb.gsfc.nasa.gov/misc/NotesW.html#WI_OR_PRE, respectively. Mars- 
Odyssey ephemerides were taken from JPL’s HORIZONS system https://ssd.jpl.nasa.gov/ 
horizons.egi 

For the near-Earth spacecraft and Wind, ephemeris uncertainties contribute less than 1 ms to the 
propagation time delay, so we conservatively assume a systematic error in Konus-GBM and Konus-BAT 
triangulations to be 1 ms. For the Konus-HEND and GBM-HEND triangulations we take 360 ms as the 
30 systematic uncertainty. 

For GRB 200415A triangulation we used the following ligtcurves: 2 ms Konus (see Konus-Wind 
section), 390-1600 keV; 0.1 ms GBM, 360-1000 keV, constructed from the TTE data of triggered detectors 
(0, 1, 2, 3, 4, 5, 9, and a; only Nal data were used); 0.1 ms BAT, 25-350 keV, constructed from the 
TTE data from the GUANO system; 250 ms HEND, 50-3000 keV; and 7.8 ms INTEGRAL-PICSIT, 
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Methods Table 1. Triangulation annuli: the instruments involved in triangulation and the lightcurve 
temporal resolution used (1st column), the annulus center right ascension and declination in the equatorial 
J2000 system (2nd and 3rd columns, respectively), the annulus radius (4th column), and its half width 
corresponding to 30 statistical cross-correlation time delay uncertainty with systematics added in 
quadrature. 


Instruments R.A. (J2000) Dec.(J2000) R OR 
involved (deg) (deg) (deg) (deg) 
GBM(0.1 ms)-KW(2 ms) 1.9406 2.0665 28.9781 0.0298 
BAT(0.1 ms)-KW(2 ms) 2.0920 2.0554 28.9243 0.0262 
KW(2 ms)-HEND(250 ms) 313.0127 -18.9203 54.5051 0.0394 
GBM(1 ms)-HEND(250 ms) 313.3351 -18.8199 54.2444 0.0391 
KW(2 ms)-PICSIT(7.8 ms) 4.1624 -2.1891 24.2751 0.1681 


250-2000 keV. 

Using these data we derived five annuli (Methods Table 1). The final IPN 30 box was constructed from 
Konus-GBM and Konus-HEND annuli (Methods Table 2). We used the Konus-GBM annulus instead of 
the narrower Konus-BAT one due to the similarity of the two instruments’ energy bands. GBM saturation 
occurred near the burst peak and does not significantly affect the cross-correlation with the Konus 2 ms 
resolution light curve. 

The annuli were combined to yield an error ellipse?^. with a major axis corresponding to the Konus- 
GBM annulus, and a minor axis corresponding to the Konus-HEND annulus. We obtain a 30 error ellipse 
centred at R.A.(J2000) = 11.885 deg, Dec.(J2000) = -25.263 deg with major and minor axes 6.25 and 
4.07 arcmin, respectively, and position angle of 61.135 deg. The ellipse area is 20 arcmin?. 

The ellipse contains the central part of the Sculptor galaxy (NGC 253; Figure 1). The image was 
obtained with the Galaxy Evolution Explorer GALEX^? during the GALEX Nearby Galaxies Survey?? 
(Observation ID: 2482083865531777024). The NGC 253 image was obtain via MAST portal (https: 


//mast.stsci.edu). 


Konus-Wind 


Konus-Wind?" (KW) consists of two identical Nal(TI) scintillation detectors, each with 27 sr field of view, 


mounted on opposite faces of the rotationally stabilized Wind spacecraft^?, such that one detector (S1) 
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Methods Table 2. The 30 IPN box: The error box area is 17 sq. arcmin, and its maximum (minimum) 
dimensions are 7 arcmin (4 arcmin). The Sun distance was ~ 37 deg. 


Box center/ R.A. (J2000) Dec.(J2000) 
vertices (deg) (deg) 
Center 11.885 (00h 47m 32s) -25.263 (-25d 15m 47s) 
1 11.846 (00h 47m 23s) -25.308 (-25d 18m 29s) 
2 11.931 (00h 47m 43s) -25.279 (-25d 16m 44s) 
3 11.923 (00h 47m 42s) -25.218 (-25d 13m 05s) 
4 11.839 (00h 47m 21s) -25.247 (-25d 14m 49s) 


points towards the south ecliptic pole, thereby observing the south ecliptic hemisphere, while the other (S2) 


observes the north ecliptic hemisphere. 


Each KW detector is a cylinder 5 inches in diameter and 3 inches in height, placed into an aluminum 
container with a beryllium entrance window. The crystal scintillator is viewed by a photomultiplier 
tube through a 20 mm thick lead glass, which provides effective detector shielding from the spacecraft's 
background in the soft spectral range. The detector effective area is ~ 80-160 cm?, depending on the 
photon energy and incident angle. The energy range of gamma-ray measurements covers the incident 


photon energy interval from 20 keV to 20 MeV. 


The instrument has two operational modes: waiting and triggered. While in the waiting mode, the 
count rates (lightcurve) are recorded in three energy band covering the ~ 20-1500 keV energy band, see 
Methods Table 3, with 2.944 s time resolution. When the count rate in G2 exceeds a ~ 90 threshold above 
the background on one of two fixed time-scales, 1 s or 140 ms, the instrument switches into the triggered 


mode. 


In the triggered mode, light curves are recorded in the same bands, starting from 0.512 s before the 
trigger time 7o with time resolution varying from 2 ms up to 256 ms. For the bursts of interest here, the 
whole time history is available with 2 ms resolution. 

Multichannel spectral measurements are carried out, starting from the trigger time 7 (no multichannel 
spectra are available before 70) in two overlapping energy intervals PHA1 and PHA2 (Methods Table 3). 
with 64 spectra being recorded for each interval over a 63-channel, pseudo-logarithmic energy scale. The 


first four spectra are measured with a fixed accumulation time of 64 ms in order to study short bursts. 
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Methods Table 3. Konus-Wind calibrations for GRB 200415A and GRB 051103 


Burst Det. Inc. angle Gl G2 G3 PHAI PHA2 

(deg) (keV) (keV) (keV) (keV) (keV) 
GRB 200415A SI 62.2 22-90 90-390 390-1600 28-1600 330—20000 
GRB 051103 S2 70.8 17-70 70-300 300-1200 20-1170 240-14800 


For this analysis we use a standard KW dead time (DT) correction procedure for light curves (with a 


DT of a few microseconds) and multichannel spectra (with a DT of ~ 42 microseconds). 


Temporal analysis For the temporal analysis we used time histories from 7 — 0.512 s to Tọ + 0.512 s 
in three energy bands: G1 , G2, and G3 with a time resolution of 2 ms . The total burst duration 7100, and 
the 799 and 759 durations (the time intervals that contain 5% to 95% and 25% to 75% of the total burst 
count fluence, respectively’), were calculated in this work using the light curve in the ~ 80-1500 keV 
energy band (G2+G3). Burst start and end times in each band were calculated at the 50 level with a 
method similar to that developed for BATSE??. The background count rates, estimated using the data from 
~ Ty — 2500 s to ~ Ty — 150 s, are 958.7 s! (G1), 349.5 s! (G2), and 223.0 s! (G3); and 1080.3 s~! 
(G1), 394.0 s^! (G2), and 135.5 s^! (G3) for GRB 200415A and GRB 051103, respectively. 


Spectral analysis For the bursts of interest we analysed both multichannel and three-channel KW energy 
spectra. The multichannel spectra accumulation intervals are presented in Extended Data Tables 2 and 3. 
The background multichannel spectra were extracted in the intervals from To + 8.448 s to To + 491.776 s 
and from Tọ + 98.560 s to To + 491.776 for GRB 200415A and GRB 051103, respectively. The emission 
evolution at a finer time scale can be explored using three-channel spectra, constructed from the counts 
in the G1, G2, and G3 energy bands in the six intervals (Extended Data Tables 2 and 3). Details on KW 
three-channel spectral analysis can be found elsewhere?! . 

We performed the spectral analysis in XSPEC, version 12.10. 1??, using the following spectral models: 
a simple power law (PL), a custom exponential cutoff power-law (CPL) parametrized with peak of v Fy 
spectrum and energy flux as the model normalization, the Band GRB function®”, a single blackbody (BB) 
function with the normalization proportional to the surface area, and a sum of CPL and BB functions 


(CPL+BB). The details of each model are presented below. 
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The power law model: 


feu = A(E/En)®; (1) 


the custom exponentially cutoff power law (CPL): 


n(E) (E/En)” exp(—E(2+ @)/Ep) 


Emax 
foL = Fxn(E)/ f n(E)EdE ; 
Emin 


the Band function: 


fea cay CIE)" ero LRS). E<(a-B) x 
Band E E B Ep(a-B) (a—p) u E > _ (Ey 
( / n) En(2+0) exp(B a), Z (a B) 2+0? 


where f is the photon spectrum, measured in photon cm? s^! keV~!, A is the model normalization, 
E, = 100 keV is the pivot energy, Ep is the peak energy of the vF, spectrum, and F is the model energy flux 
in the Emin — Emax energy band; œ and p are the low-energy and high-energy photon indices, respectively. 


The single blackbody (BB) function is the bbodyrad XSPEC model. 


The Poisson data with Gaussian background statistic (PG-stat) was used in the model fitting process as 
a figure of merit to be minimized. The spectral channels were grouped to have a minimum of one count 
per channel to ensure the validity of the fit statistic. Since the CPL fit to a three-channel spectrum has 
zero degrees of freedom (and, in the case of convergence, PG-stat=0), we do not report the statistic for 
such fits. The 68% confidence intervals of the parameters were calculated using the command steppar 


in XSPEC. 


A summary of constrained spectral fits with CPL, BB, and CPL+BB models is presented in Extended 
Data Tables 2 and 3. For both GRB 200415A and GRB 051103, the PL model failed to describe the spectra, 
with PG-stat/dof >10 in all cases. Use of the Band GRB function does not constrain the high-energy 
photon index D for GRB 2004154 spectra, and only marginally improves the CPL fit to the time-integrated 


spectrum of GRB 051103, with similar, within errors, Ep and a, and po-3. 
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Burst energetics For both bursts, the total energy fluence S was derived using the 20 keV-10 MeV 
energy flux of the best-fit (CPL+BB) spectral model. Since the time-integrated spectrum accumulation 
interval differs from the 7100 interval, a correction which accounts for the emission outside the time- 
integrated spectrum was introduced when calculating S. 

The peak flux Fpeak was calculated on the 4 ms scale using the energy flux of the best fit with CPL 
model to the three-channel spectrum at the peak count rate interval (Tọ — 0.002 s-70 + 0.002 s). We note, 
that the peak flux of GRB 051103 estimated in this work is a factor of ~ 2.5 lower than that reported 
from the previous analyses of KW and RHESSI data”, which used wider spectral intervals and did not 
separate the relatively soft spectrum in the huge 4 ms spike (Ep ~ 1.2 MeV) and the considerably harder 


emission observed immediately after its falling edge (Ep ~ 3 MeV). 
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Data Availability 


Links to the Wind ephemeris and clock accuracy data are provided in the Methods section. The Konus-Wind 
lightcurve and spectral data are available at Ioffe web site http://www.ioffe.ru/LEA/papers/ 


SvinkinNat2020/data/. The HEND lightcurve is available from the HEND team on reasonable 


request. The Fermi, Swift, and INTEGRAL data are freely available on-line. 


Code availability 


XSPEC is freely available on-line. 
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Extended Data Figure 1. GRB 051103 (red stars) and GRB 200415A (blue stars) as possible 
cosmological GRBs at different redshifts (0.01 < z < 1). The Konus-Wind samples of short/hard and long 
GRBs with known redshifts?? are shown by triangles and circles, respectively. The recent update?? for the 
hardness-intensity relation in the cosmological rest frame (Ep z—Ejso, ‘Amati’ relation) is plotted with the 
solid line together with its 68% and 90% prediction intervals (dashed lines). Considering only its 
spectrum and energy fluence, GRB 200415A is consistent with the KW sample of short GRBs if at 
redshift z ~ 0.05 — 1. In the case of GRB 051103, the implied short GRB redshift is z ~ 1, and intrinsic 
Ep ~ 5 MeV. 
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Extended Data Table 1. Summary of GRB 200415A and GRB 051103 properties 


Parameter GRB 200415A GRB 051103 


Host galaxy (distance) NGC 253 (3.5 Mpc)  M81/M82 group (3.6 Mpc) 


Temporal properties 


Trise (ms) <2 <4 

Ter (ms) ~ 50 ~ 50 
tgux (ms) ~ 30 ~ 30 
Too (8) 0.138 0.324 

Too (s) 0.1002: 0.014 0.138 + 0.020 
Tso (s) 0.048 + 0.005 0.058 + 0.004 


Peak spectrum 7p (-0.002 s, 40.002 s), CPL model 


CPL photon index 0 —0.59*2.17 —0.13*015 
CPL Peak energy E, (keV) 1190+460 1250+350 


Time-integrated spectrum Tọ (0, +0.192 s), CPL+BB model 


CPL photon index & —0.0210.38 0.081) 78 
CPL Peak energy Ep (keV) 1080+710 2690+210 
Blackbody temperature kT (keV) 99+37 107414 
Blackbody radius R (km) 23339 (@3.5 Mpc) 37*6 (@3.6 Mpc) 
Blackbody contribution to flux ~14% ~9% 


Peak energy fluxes (erg cm7? s_!), in the 20 keV-10 MeV band 


4 ms scale, To (-0.002 s, +0.002 s) Oe rex 1074 11375 x 1074 
16 ms scale, To (-0.002 s, 40.014 s) Lir x 1074 8.908757 x 1074 
64 ms scale, Tọ (-0.002 s, 40.062 s) 0.43*007 x 10-4 4,381 76 x 1074 


Energy fluences (erg cm ?), in the 20 keV-10 MeV band 
Initial spike 3886 o x 10 5 4.61139 x 1076 


To (-0.002 s, +0.002 s) To (-0.002 s, +0.002 s) 


Total Ra x 1076 34.3122 x 1076 
To (-0.004 s, +0.192 s) Ta (-0.006 s, +0.192 s) 


Flare entergetics, in the 20 keV-10 MeV band 


Liso, 4 ms scale (erg s^!) ^. 1.4 x 1088 ~ 1.8 x 1048 
Eiso (erg) ~ 1.3 x 1046 ~ 5.3 x 1046 
KW maximal detection distance (Mpc) ~ 13.5 ~ 15.8 
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Extended Data Table 2. GRB 200415A spectral fits 


Time interval (s) model a E, kT R Flux (20 keV-10 MeV)  PGstat/dof 
(keV) (keV) (km) (1076 erg cm"? s7!) 


Three-channel spectra 


-0.004 — -0.002 CPL —028' 5 5201320 - - sit - 
-0.002 — 0.002 CPL = -0.59'07 11907200 — — = 9601730 - 
0.002 — 0.032 CPL 0307075 980712 - - 111471 - 
0.032 — 0.064 CPL 0.164922 7104160 - - 437 - 
0.064 — 0.096 CPL 0441090 — 496194 = - 19.3777 = 
0.096 — 0.128 BB - - 795. ad 6.81) 0.5/1 


Multichannel spectra 


0.000 — 0.064 CPL 012: 5- “106% = = 85.4163 55/64 
0.064 — 0.128 CPL 030707. 458178 - - 12,5725 49/47 
0.128 — 0.192 BB - - 71 265, i405 22/31 
0.000 — 0.192 CPL 0.014912 ay = = 32.3+24 67/75 

CPL+BB  —0.02'023 10801710 99t3} 23*16 33.3121 63/73 
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Extended Data Table 3. GRB 051103 spectral fits 


Time interval (s) model a Ep kT R Flux (20 keV-10 MeV)  PGstat/dof 
(keV) (keV ) (km) (1076 erg cm? s7!) 


Three-channel spectra 


-0.004 — -0.002 CPL | —032'02 1380289 - = 2074817 - 
-0.002 — 0.002 CPL Od 1250+500 — n - 1150+540 - 
0.002 — 0.032 CPL 020: 19.  3620+7280 - - 9407128 - 
0.032 — 0.064 CPL 0.644973 — 9307250 = = 11675 - 
0.064 — 0.096 CPL 02802 — ort - - 42*7 - 
0.096 — 0.128 BB - - 935 Tu 13.1+]$ 0.8/1 


Multichannel spectra 


0.000 — 0.064 CPL — —dup os. 257015 = = 444437 93/77 
CPL+BB 03010) 299000. 120" 35:505) 444432 83/75 

0.064 — 0.128 CPL GATTI 500° - - 34.2124 57/58 
0.128 — 0.192 CPL 0i «320% - - 5.5419 29/43 
BB = = 789 «Asl Si 30/44 

0.000 — 0.192 CPL = —-0.30'000 2300+150 - - 16279 127/85 
CPL+BB 0.087} 75 2690+70 + 10771, 36° 16271) 98/83 
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